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Abstract 

We extend two of the methods previously introduced to find discrete 
symmetries of differential equations to the case of difference and differential- 
difference equations. As an example of the application of the methods, we 
construct the discrete symmetries of the discrete Painleve I equation and of 
the Toda lattice equation. 



1 Introduction 

Symmetries have always played a very important role in the study of differential 
equations Lie symmetries have been introduced by Sophus Lie as a tool to 
unify all solution techniques for ordinary differential equations. In particular they 
are useful to 

• get special solutions by symmetry reduction. 

• classify equations according to their symmetry. 

• prove equivalence of equations under point transformations. 
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They have been extended with success to the case of differential-difference 
and difference-difference equations. As was shown in [2| the invariance of a 
differential-difference equation with respect to a shift of the lattice h = n + N 
provide a reduction of the equation to a system of N ordinary differential equa- 
tions. Moreover, discrete symmetries, i.e. symmetries associated to a discrete 
finite group, are very important in Quantum Mechanics. In this field one speaks 
of parity, charge conjugation, rotations by n, etc. . . and one uses discrete symme- 
tries to provide selection rules. 

Discrete symmetries are usually easy to guess but difficult to find in a system- 
atic way. They can be obtained by considering the normalizer of the continuous 
Lie point symmetries of the equation. This is a well known technique and can 
be found in many textbooks (see for instance O). The application to the case of 
differential equations has been considered in all details by Hydon |7|. The nor- 
malizer can be obtained only in the case when a nontrivial non commuting group 
of Lie point symmetries exists. 

In 1996 Gaeta and Rodriguez ||8|| introduced a modification of Lie technique 
which allows us to find discrete symmetries even when the Lie group is trivial. 
This technique, as we will see in Section 3, depends on a starting ansatz and thus 
it does not provide a complete result. 

Reid, Weih and Wittkopf [3| solved the determining equations for the group 
transformations directly to get the discrete symmetries. In more recent works II4I5I1 
they simplified the work taking advantage of the infinitesimal results; however, the 
method requires extremely heavy calculations. 

In Section 2 we briefly summarize the steps necessary to obtain Lie symme- 
tries for difference equations while in Section 3 we present two methods intro- 
duced to obtain discrete symmetries. As is evident from Section 2 all the steps 
necessary to obtain the discrete Lie symmetries are not significantly different in 
the case of equations on the lattice from the continuous case. The only essential 
difference is in the prolongation formula. So in Section 4 we just present two 
examples of the derivation of discrete Lie symmetries for discrete equations. In 
Section 5 one can find some concluding remarks. 

2 Point symmetries of discrete equations 

Let's briefly summarize the steps necessary to obtain Lie point symmetries for 
difference equations. More details can be found in [9l'T0'|. 

We will consider just the case of a scalar difference equation in two indepen- 
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dent lattice variables as this case will cover all the examples considered in Section 
4. The case of a differential-difference equation is obtained from the difference- 
difference case by carrying out the continuum limit in one variable. 

(i) 

A discrete equation of order Ni in a discrete variable Xn'm, i = 1,2, is a 
functional relation between Ni + 1 points in the lattice of the variable Xn^m- To 
be able to solve the discrete equation apart from the functional relation we need 
to know how the lattice points are defined. This implies that we need four further 
equations in the case of two independent variables and one for one independent 
lattice variable. Some of these equations may be trivial if the lattice is orthogonal 
and with constant spacing, but still are necessary to fix the symmetries and the 
solutions of the difference system. 

So we have: 

^{•^n+j.m+iy '^n+j,m+i) = 0, (2.1) 
^a.i.-^n+j,m+ii'^n+3,m+i) = 0, (2.2) 

Un,m = u{x^n^^, x^^^), 1 < a < 4; k = l,2; -ii < i < 12] 

-ji <3< 32] «i,«2, ji, j2, e z-°, 

where A is the difference equation, Ea are the equations determining the two 
independent lattice variables and n and m are some indices which characterize 
them. System d2. 112.21) must satisfy some obvious conditions so as to be able to 
calculate the variables in all the lattice plane. 

A Lie point symmetry for equations (12.112.21) is a point transformation 

'^n,m ~ -^A i.'^\i]mi'^'ri,m)i '^n,m — ^ \{-^n)ni '^'n.,m) j ~ l)^), (2.3) 

which leaves equation ( 12.112.21 ) invariant. The symbol A indicate the group param- 
eters. As this transformation acts on the entire space where the independent and 
dependent variables are defined, the functions F^'' and G\ determine the transfor- 
mation everywhere. We can introduce the corresponding infinitesimal transforma- 
tions of coefficients ^^"^^xH'Jn, Un,m), 4>{xn,L, Un,m) and thus write the vector field 

-^n,m = ^^\xl^lny'^n,rn)d H) + 0(2^1 Mn,m)f^M„m ; (2.4) 

and its prolongation 

32 12 

Pr^n,m = ^ ^ ^ ^ ^n+l,m+j- (2.5) 

j=-ji l=-ii 
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The invariance conditions, i.e., the necessary conditions which provide the 
symmetries for equations ( 12.112.21 ). are given by 

prX„,^A = when (A, = (0, 0), (2.6) 



WXn,rnEa = Q whcn ( A, E„) = (0, 0) . (2.7) 

Equations (12.6112.71) are a set of equations for and 0. In equations (12.6112.71) 
and <p will appear in various points of the lattice, which, when the equations 
A = and Ea = Q have been taken into account, are all independent. So we have 
five functional equations for the functions and cj). The variables appearing 
in these five equations are all independent. These independent variables appear 
either explicitly in the equation or in the unknown infinitesimal coefficients. As 
the infinitesimal coefficients are analytic functions we can convert the determin- 
ing equations into a system of differential equations by differentiating them with 
respect to the independent variables. In such a way we get an overdetermined sys- 
tem of, in general, nonlinear partial differential equations. We solve the obtained 
equations and introduce its solution into the functional equations and solve them. 

As an example of this procedure we will calculate the continuous symmetries 
of a discrete Painleve I equation lITTIl . The discrete Painleve I equation is given by 

axn + /3 

Un+l +Un + Un-1 = h 7 (2-8) 

Un 

where a, (3 and 7 are arbitrary constants and m„ = The equation which 

defines the lattice points can be written as: 

Xn+l - Xn = h (2.9) 

so that Xn = hn + xq. With the choice Un = —I + h'^w{x), a = —h^, (3 = 3 and 
7 = 6 equation (12.81) reduces, in the continuous limit, to 

w^, = 6w^ + x. (2.10) 



which is the Painleve I transcendent [T2||. 

The infinitesimal symmetry generator is given by 

Xn = in{Xn,Un)d^,, + , , (2.11) 

and its prolongation 

1 

WXn,m = Xn+j- (2.12) 

i=-i 
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Applying eq. (12.121 ) to the lattice equation ( 12.91) we get 

^n+l (-^n+l ) ^n+l) ^n(^-^ny ^n) : 



(2.13) 



which, taking into account that m„ and Un+i are independent variables, gives that 
cannot depend from u„ and, as a function of n, must be a constant (cannot vary 
between different points of the lattice) 



(2.14) 



In a similar way, by applying (12.121) to equation (12.81) and differentiating the result 
with respect to m„ and Un+i, one can prove that 



(2.15) 



where 0n and 0n "* satisfy the following overdetermined system of equations 



(aa;„ + /3)0(f) 



/,(0) 



+ 



r'n ' 

aKo, 

0, 

0. 



(2.16) 
(2.17) 
(2.18) 
(2.19) 



According to the values of the parameters (a, P and 7) we have various possi- 
bilities: 



If « 7^ then from equation ( 12.171) Kq = and 



(Em 



0. From equation 



and thus no symmetry is present. 



If a = but /3 7^ then from equation (12.171) . (pn^ = and from equation 
(irT8l) . 0^°^ = 0; so only Kq ^ 0, i.e., only space translations are possible. 

If a = 0, /9 = equations (12.171 12.181) are identically satisfied, equation 
(irT6l) implies that (pll^ = Ki and ifTSll 





n 


-^(iV3 + l) 









(2.20) 



So in this case the now linear equation has a four dimensional symmetry 
group. 
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To end this section let us consider the case of symmetries of differential- 
difference equations. The passage from a difference-difference equation to a 
differential-difference equation is done by carrying out the continuum limit for 
a lattice variable, say Xn,L = tn,m, when the distance between the points along 
this direction goes to zero and the lattice index goes to infinity in such a way that 
the position t remains finite. This implies that the corresponding lattice spacing 
r cannot be modified by a point transformation but it is a fixed number which 
can tend to zero. So, for example, the lattice variable tn,m cannot be described by 
dilation invariant equations like 

tn+l,m '^tn,m ~l~ ^ra— l,m 0, tj^ m+1 tn,m 0, (2.21) 

as in this case the parameter r is an integration constant and not a parameter of 
the equations. Let us choose, consequently, the lattice equations for the variable 

tn,m 

^n+l,m ^n,m "^i ^n,m+l ^n,m 0. (2.22) 

The solution of equation (12.221) reads 

t = tn,m = Tn + to, (2.23) 

where, in all generality, we can set to = 0. As for the remaining lattice variable 

(2) 

Xn'm = Xn,m, its position is not changing in time, i.e. 

Xn+lj-m 2^71, m 0, (2.24) 

while its variation along the lattice, i.e., along m, can depend on t and on Un,m in 
a way that is defined by one of the equations (12.21 ). say Ei = 0. 

The continuous limit is obtained by considering n ^ oo as r — 0. In such 
a limit the variables Xn,m and Un,m will no more depend on n, equation (12. It . 
if it has the proper dependence from r, will reduce to a differential-difference 
equation for Um = u{t, Xm) and equations (12.21) will reduce to just one equation 
for the lattice variable Xm, Ei = 0, while the other equations (12.221 12.241 ) are 
identically satisfied in the limit. 

In this limit the symmetry vector X„ given by equation (12.51 ). reduces to 
Xm, which will inherit the properties of Xn,m as applied to the lattice equations 
(I2.22ll2!24l) even if these equations in the limit reduce to = 0. Applying equation 
(12.61 ) to equations (12.221 12.241 ). with Xnjn = tn,m and Xnln = Xn,m, we get the 



6 



following three determining equations 





(2.25) 
(2.26) 
(2.27) 



As the differential-difference equation will involve at least m^+i Un,m., Un.m+i, 
we can always take Un.m, Un,m+i as independent variables and express Un+i,m in 
term of them. By differentiating equation (12.251) with respect to u„ m+i we get 
■C^^n^n,™, Xn,m)- By a similar reasoning for the variable Xn,m we can reduce the 
function ^^^^ to ^^^\tn^m), the function ^^^^ is just a function of t. In a similar 
way we will find that 

and thus, in the continuum limit, we must have 



X„, = ^^'\t)dt + ^^^\t,Xm)d,^+(p{t,x^,Um)du^. (2.29) 

3 The calculus of discrete symmetries 



We will now discuss two methods of determining the discrete symmetries of dif- 
ferential equations. One of them was proposed by Hydon QIHi and it is essen- 
tially the classical method of constructing the normalizer of a group. The other, 
due to Gaeta and Rodriguez |8|, is a modification of Lie's method, defining a 
discrete symmetry as a discretization of the parameter of a continuous symme- 
try. There are at least two other methods to construct discrete symmetries. One 
of them is based in a formulation of differential equations through differential 
forms fT5\ and the other [3,4,5| solves the determining equations for the group 
transformations directly. They will not be considered in this work. 

3.1 Automorphisms of the symmetry algebra 

Let us consider a differential equation for a dependent variable u and N inde- 
pendent variables x = (xi, . . . , xn), with a Lie group of symmetries, G and its 
corresponding Lie algebra Q. We construct all the automorphisms of this Lie alge- 
bra. Some of them will, obviously, correspond to continuous symmetries. Others 
will be essentially new and will define discrete symmetries of the equation. And, 
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finally, others will not be symmetries of the equation. These ideas are very well 
known in the theory of Lie algebras and groups. Inner automorphisms correspond 
to conjugation by elements of the group, while outer automorphisms have not this 
character. See for instance [6| where the complete and general Lorentz group are 
obtained in this way from the proper Lorentz group. 

Let Qhe a Lie algebra of finite dimension n and B = {Xi, . . . , X„} a basis of 
Q. The commutation relations in this basis can be written as: 

[X„Xj]=c1^Xk, (3.1) 

where cf^ are the structure constants of Q (a sum is understood over repeated 
indices running from 1 to n). 

The defining equation for an automorphism (p : Q — > Q is 

0[X,r] = [0(X),0(F)], (3.2) 

which, in the basis B, reads: 

= 40(Xfc) = [0(X,),0(r,)]- (3.3) 
Let $ be the matrix representation of (j) in the basis B (det $ ^ 0). Then 

= ^\X^ (3.4) 

and eq. (13.31) is written as: 

'^ij^ k ^ j^mr (3-5) 

for all indices I (the equation is skewsymmetric in If we define the ma- 
trices of the adjoint representation 

C{^r, = 4, (3.6) 

equation (13.51) can also be written as: 

$C(z) = $'.C(/)$. (3.7) 

Automorphisms which are related through conjugation by an element of the sym- 
metry group will be considered equivalent. Thus the matrix $ can be simplified 
by conjugations with the symmetry transformations (at least if the algebra is not 
Abelian). This can be done in terms of the adjoint representation. 
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Let us consider those elements gi G G which are generated by an element 
of the basis B of the corresponding Lie algebra Q, 

^^ = e"^^ (3.8) 

The conjugation provides the transformation 

(f)-^e-^^^(f)e^^\ (3.9) 

We can consider the transformation of the elements of the basis B and we have: 

e-^^^Xj e^^^ = A{i, A)^Xfe. (3.10) 

It is not difficult to show that 

A(i,A)=e^^«. (3.11) 
We can apply the transformations (13.1111 : 

<|,^e^^(^)<l> (3.12) 

to simplify the matrix $. 

We can now construct a representation of the automorphism in the space of 
variables. Let X be a vector field given by (the indices a, b run from 1 to N) 

X = r(x,n)9,, +<^(x,n)9„ (3.13) 

which is a symmetry of the equation under study, and let us consider a symmetry, 
given by the transformation: 

(x, u) — * (x(x, m), ^(a;, u)) (3.14) 

The vector field X is transformed into a new vector field 

X = e{±,u)d,^ + 0{ic,u)di,. (3.15) 

As the transformation is a symmetry of the equation, X must have the same form 
in the new variables: 

X = C{^,u)d,^ + ^{^,u)d^. (3.16) 
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If we consider a basis of the symmetry algebra, {Xj}, i = 1, . . . ,n, the trans- 
formed vector fields are: 

Xi = ^^{±,u)d,^+iPi{±,u)d^ (3.17) 
and, as the transformation is an automorphism of the algebra, 

X, = ^^,Xj. (3.18) 
Applying (13.181) to the new variables, we get: 

Xda = (X, U) = ^\X,Xa = (e)(x, u)^^+ (x, u)^^ (3.19) 



and 



X,u = $,(x, u) = <l>\X^u = <l>^, (^el(x, + ¥^,(x, «)|^ j (3.20) 



or: 



e)(x,«)^ + ^,(x,«)^ = (3.21) 
el(x,«)^ + ^.(x,n)|^ = (*-^)Xx,^). 

We have to solve eqs. (13.211) to find the expression of the automorphism as a 
transformation in our space of variables and functions. After, we must check if 
the automorphism is a symmetry of the equation and it does not correspond to a 
continuous symmetry. 



3.2 Determining equation for a discrete symmetry 

We will now consider the method of the discretization of the parameter of a con- 
tinuous transformation. Let us briefly review how the method works. To simplify 
the description (see [iSl for a detailed exposition), let us consider a smooth curve 
in R^, y = f(x) and a vector field 

X = ax,y)d, + ^{x,y)dy (3.22) 

The point {x,y) is transformed under the infinitesimal action as 

x' = X + X^{x,y), y' = y + X(p{x,y) (3.23) 
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and the graph of our curve is transformed into a new one, defined by the trans- 
formed function: 

/a(x) = fix) + \Mx,y) - ax,y)dj{x)] (3.24) 

If we introduce a function F(x] A) such that F(x; A) = fx{x), F, taking into 
account (13.241) will satisfy the partial differential equation 

+ ^(a:, F{x; ^))^^^ = ^(^, F{x; A)). (3.25) 

In order to recover the original function, we have to impose the initial condition: 

F{x-0) = f{x) (3.26) 

and if we want to obtain the same function for a particular value Aq of the param- 
eter and hence a discrete symmetry 

F(x;Ao) =F(x;0) = /(x) (3.27) 

This is equivalent to find periodic solutions of the equation (13.251) . This equation 
is a functional equation and hard to solve. 

If we consider a differential equation instead of the graph of a function, we 
have to pose the same question in the appropriated jet space [ 1 1. The answer is 
simpler in this case, as the right hand side of equation (13.251) is computed from the 
prolongation of the vector field under consideration. 

Let us consider a differential equation: 



djU = /(x,M,(9j'M, . . .) 



(3.28) 



where dju = di\ . . . di^u, J = (ji, ...Jn)- 

The determining equation, used to describe the discrete symmetries, is then: 



^ ^ dxi ^ dujt 



dX 



(3.29) 



dju=F 



where (f)"^' are the J'-prolongations of the vector field (13.131) and one has to rewrite 
equation (13.281) in terms of F. Then one looks for periodic solutions of this equa- 
tion in A. 

For instance, consider a differential equation of the following type: 



Utt = f{^,U,Ut). 



(3.30) 
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Choosing F 



F(x, u, Ut, A), the determining equation is: 




dttu=F 



(3.31) 



and 



F(x, u, Ut] 0) = F(x, u, Ut] Ao) 



f{^,U,Ut) 



(3.32) 



for some Aq. 

To apply this method we do not need to know the continuous symmetries. In 
fact, the determining equation provides a general solution which does not depend 
on the differential equation we are studying. The equation appears when we im- 
pose the boundary condition as in (13.261 ) and the discrete symmetries by requiring 
equation (13.271) . However, we have to make some assumptions to get a solution of 
the equation, because we do not know the symmetries which appear in it. We will 
use this method to determine the discrete symmetries of the Painleve I equation 



4 Lie discrete symmetries of lattice equations 

The construction shown in Section 2 of point symmetries of discrete equations 
is very similar to the standard approach to point symmetries of continuous equa- 
tions [il. The main difference lies in the form of the prolongation and in the 
way we solve the determining equations. Consequently, the procedure presented 
in Section 3 for constructing discrete symmetries of continuous equations can be 
carried over in a straightforward manner to the discrete case, just by changing the 
form of the prolongation. So, in the following we will just apply the methods 
discussed in Section 3 to two examples of equations of Mathematical Physics, the 
discrete Painleve equation (12.81) and the Toda equation 1 16|. The Painleve equa- 
tion, as we have shown in Section 2, has no continuous symmetry (for a generic 
choice of the parameters), and we will use the determining equation method to 
find its discrete symmetries. In the Toda equation we have a continuous group of 
symmetries II17LI18I and we will construct the discrete symmetries by computing 
the automorphisms of the corresponding Lie algebra. 

4.1 Discrete Painleve equation 

As we have said in the introduction to this section, the discrete Painleve I equation 
( 12.81 ) has no continuous symmetry so that we cannot apply the method of the 
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normalizer. 

We define the extended A dependent equation in the transformed space as 

u+ + u_ = F{x,u;X), (4.1) 

where, for A = 

F(x,m;0) = + + (4.2) 

u 

For the sake of simplicity we have used the notation 

U = u{Xn), U+ = U- = u{Xn-l)- 

The generic Lie point transformation, written in terms of the infinitesimal sym- 
metry generators, is: 

dx du , / ^ ^\ 

— = ^{x,u), —=(f){x,u). (4.3) 



We consider the same definition of the lattice (12.91 ) that we used for constructing 
the Lie point symmetries in Section 2. So we have ^ = Kq. Differentiating 
equation (14.11 ) with respect to A and taking into account equation (14.31 ) we get 
(ED: 

dF dF OF 

— + Ko— + (j){x,u)— = (j){x+,u+) + (j){x_,F-u+) (4.4) 



By differentiating equation (14.41) with respect to and u we get 

(f){x,u) =(j)^'^\x)+u(j)^^\x). (4.5) 

Substituting equation (14.51) into equation ( 14.41) and requiring that the obtained 
equation be satisfied for any m+, we get 

= ^1 + K^i-l)^. (4.6) 

Equation (14.41) is now reduced to a linear partial differential equation of first order 
which can be solved on the characteristics: 

dA dx du dF 
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where 



(0 



0,1. 



(4.8) 



The first invariant is: 

x = x + KoX. (4.9) 

The action of this transformation on the discrete Painleve equation (12.81) corre- 
sponds just to a change in the parameter p. We could carry out the calculation 
with Kq ^ 0, but for the sake of clarity of the presentation, we have set i^o = as 
this transformation cannot provide any discrete Lie transformation. Consequently, 
will not depend on A. 

The next invariant is obtained by integrating equation (14.71) for n as a function 
of A. We have: 



u 



U + 



1-e 



(4.10) 



The integration of equation (14.71) for F satisfying the boundary conditions (14.21) 
gives 



F{x, u; A) 



(1) 



0(1) 



:i-e 



+ 



ax + P 



+ 7 



(4.11) 



To get a discrete Lie symmetry we require that there exists a value of A, say 
Ao, such that 

F{x,u;0) = F{x,u;Xo). (4.12) 
If we want (14.121) to be satisfied, we need 







(4.13) 
(4.14) 



From (14.131) we obtain K2 = and then, (p^^^ = Ki. Equation (14.141) is solved 
by requiring one or the other of the following two conditions: 



XoKi = 27riiV, N eZ 
cj)^'^ = 



(4.15) 
(4.16) 
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In the case of condition ( I4.15I ). equation (14.121) is satisfied. However, this provides 
no discrete symmetry. 

Let us go over to the second condition (14.161) . If 7 7^ 0, equation (14.121) 
implies equation (I4.15I) . i.e., no discrete symmetry is present. If 7 = 0, then we 
have KiXq = inN, providing a discrete symmetry u = ±u even in the case when 
no continuous symmetry is present. However, in this case, the continuum limit of 
this difference equation is not Painleve I. 



4.2 Discrete symmetries of the Toda equation 

Let us consider as our second example the Toda equation: 



Utt 



(4.17) 



where x± = x ± h, u± = u{x ± /i, t) and h is the lattice step (see equation (12.91) . 

As it is well known lfT7lfT9lfT8l . the following operators form a basis of the 
symmetry algebra of the Toda equation: 

2r 

Xi = 9„, X2 = 5,, X^ = du X, = tdu, X,=tdt-—d^ (4.18) 

h 



The nonzero commutation relations are: 



— Xi, [X3,X4] — Xi, [XsjXs] — X3, 



[X2.X,]- 

and the matrices C{i) of the adjoint representation are: 



-X4 
(4.19) 



= 0, C(2) 



-2/h \ 



C(3) 



1 ■ \ 
■ 1 



2/h 



C(4) 



C{5) 



(4.20) 



-1 



V 



V 



Applying equation (13.71 ) we get : 

hi = hs = hi = hi = hi = hi = ha = hi = 0, 



(4.21) 
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where, to simplify the notation we have written 
remaining equations are: 



&34(&55 + 1) = &44(&55 

bu = bssb^ - bs^b^s- 



bij. The most significant 

1) = &43(&55 + 1) = 0, (4.22) 

(4.23) 



The determinant of $ must be different from zero (0 is an automorphism) and 
consequently, taking into account equation (14.211) ). we must have bu ^ 0. Using 
equation ( 14.221 ) we conclude that 655 = ±1. We will distinguish two cases 

a) &55 = 1 In this case, 634 = 643 = and 633 7^ 0, 644 7^ 0. The other elements 
of the matrix $ must satisfy the equations: 



&32 

bl3 

bu 

The matrix $1 = ^{b^^ = 1) is: 



$1 



b) 655 = -1 Now, 633 = 644 
satisfy the equations: 



^42 = &52 = 0, 
^33^45, &14 = &3,^&. 
^22 = ^33 ^44- 



35 "44; 



(4.24) 



^ ^33^44 


bi2 


&33&45 


^35 ^44 


bi5 


\ 





^33^44 








&25 










^33 





&35 













644 


&45 




V 











1 


/ 



(4.25) 



and 634 7^ 0, 643 7^ 0. The other elements 



(4.26) 







^32 


= &42 = &52 = 0, 












bl3 


= ^^35 ^43) ^14 = 


= -6346 


45; 








bu 


= — &22 = — ^34^43 








The matrix $2 






-1) is: 










/ 


— 634^43 


bl2 ~^35^43 ~ 


634645 


615 


\ 









634643 





625 




$2 = 










634 


635 











643 





645 






V 











-1 


/ 



(4.27) 
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To simplify the matrix $ we will conjugate the automorphism using the continu- 
ous transformations in the adjoint representation. 

The exponentials of the matrices C{i) are easy to find. Using C(3) we can put 
&35 = 0, with C(4), 645 = 0, C(5) gives 612 = 0, and, finally using C(2), 615 = 0. 
Then, the simplified $1 is: 



(4.28) 





/ b 


33^44 





\ 













^33644 















633 







b33, &44 7^ 









644 











\ 








1 ) 






The same procedure can be used with $ 


2 and the result is: 




( - 


^'34^'43 








\ 











634^43 















634 







, ^'34, ^'43 7^ 









643 











\ 








-1 


) 





(4.29) 



The following step is to realize the automorphisms in the space of the variables 
and functions of the Toda equation. We have to solve the following system of 
equations: 



^j{t:X:'^)-Q^ + ^j{t:X,U)— + (pj{t,X,U) — 

dx Ox Ox 

Tj {t,X,u)—+ {t,X,u)— + ipj it, X,U) — 

Ou on on 

^3i't^X^u)—+^j{t,X,u)— + ^j{t,X,u) — 



{^-yjTi{t,X,u), 

(<^>-l)^e.(^,£,«), 
{<^-^y.(fii{i,x,u), 



(4.30) 



where j 



, 5. In Case a). 



/ 













\ 







III/ 








a 





























u 





V 














1 / 



Then, for j — 1, Xi — du 



dt ^ dx ^ du 
du ^ du ' du 



/X, i/ 7^ 0. 



(4.31) 



(4.32) 
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Forj = 2,X2 = a 



X 



For j = 3, X3 = dt 
For j = 4, X4 = tdu 



di dx du 

— = 0, T^ = ^ly. Tr = 0- (4-33) 

ox ox ox 



dt dx ^ du ^ , . ^ 



dt ^ dx ^ du ^ 
(7M a-u du 



FoTj = 5,X, = tdt-fd^ 

^di 2x di ^ ^dx 2x dx ^du 2x du 2x (4 36) 

dt h du ' dt h du ' dt h du h 
The solution of equations (14.321) . (I4.33I) . and (I4.34I) is; 

i = fit + a, X = /iz/x + /5, u = yuz/u + 7. (4.37) 

Substituting in ( 14.351) we get a = and in ( 14.361 ) we get a = and P = 0. The 
translation in m is a continuous symmetry so the possible discrete symmetries are: 

i = cit, X = C2X, u = C2U, ci, C2 7^ 0. (4.38) 

Finally, we check if this transformation is a symmetry of the equation. To do 
so, we apply it to the Toda equation (14.171) . The Toda equation in the new variables 
reads: 

—Uff = 6^=2 L V ' <=2; ^ ' _ 6^=2 L V ' ; ^ ' C2^j (4.39) 
C2 

with h = c^h. Then, C2 = ±1, ci = ±1 are the only admissible solutions. When 
C2 = 1, ci = —1, we obtain a discrete symmetry, given by the transformation: 

i = -t, X = X, u = u. (4.40) 

When C2 = — 1, ci = 1, we get another discrete symmetry: 

i = t, X = —X, u = —u. (4.41) 

It is easy to check that the automorphism of Case b) cannot be realized in the 
representation under consideration. 
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5 Conclusions 



In this article we have shown that the two methods, the automorphisms of the 
symmetry algebra 1 14| and the determining equation for discrete symmetries [(S||, 
provide discrete symmetries, even in the case of discrete equations. It is worth- 
while to notice that in the case of the Toda lattice, we have obtained a non obvious 
discrete symmetry using these techniques. 

There are some drawbacks for both methods. In the case of the automorphism 
method, if the symmetry group is very large, the matrices involved become big 
and the final defining equations are very overdetermined and, in some cases may 
need symbolic manipulation programs to carry out the calculations. Moreover, the 
method is not applicable if there are no continuous symmetries, as it is the case, 
for example, of the Painleve equations. 

In the case of the determining equation method, the equation for the function 
F can some times be undetermined. For example, in the case of the Volterra 
equation 

ut = u{u+-u^), (5.1) 
having chosen F = F{t, u, u+, A) = Ut the determining equation to solve is 

OF dF dF OF dF 

ax + ^("ar + ^<*- "'a: + ^C-"- W ^ "-'ST = + " ^'l^- 

(5.2) 

We have no hint of the form of 0(t, u) and, consequently, equation (15.21) is not 
solvable. Of course, one could try to introduce some ansatzs for this function and 
maybe find some discrete symmetries. However, as one can easily show using 
the automorphism method, Volterra equation has no discrete symmetries and an 
ansatz will not provide any conclusion. 

Work is in progress to apply the other techniques mentioned at the beginning 
of Section 3 and for the construction of lattices and discrete equations with pre- 
scribed discrete symmetries, a problem of interest in Chemistry and Quantum 
Mechanics. Another topic which attracts our attention is the study of solutions of 
discrete equations invariant under a discrete symmetry group. 
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